ARKIV FOR FYSIK Band 18 nr 1 


Communicated 13 April 1960 by O. Kirn and E. RupBERG 


On the relativistic photoeffect in the K-shell 


By Benet C. H. Nace 


With 7 figures in the text 


CONTENTS 
Page 
LEVERS OE, 2 A eR Ss Ae Se 1 
1. The expression for the differential cross section, including polarizations 3 
2. Initial and final state wave functions........5.+.+...e4e8-. 6 
3. Expansions of the reaction amplitudes in radial matrix elements . . 8 
4. An approximate final state wave function in a Coulomb field... . 10 
5. The reaction amplitudes to relative orderaZ........... #18 
Pe corbin ne IbGuSsOGarane tem: iste. ety th iS) os se 18 
a. Angular distribution, unpolarized photons............ I 
b. Angular distribution, linearly polarized photons. ........ 16 
Grilectron POlarizatvion . a) 6 sje bo sus Per re Uh 
COTE IECTORSBeCCIONU me tere Pes piers « - e le t e w LB 
Peiomianeh-ovleroyelimiluns se 2 ee es ew ee ew tl  D 
H (leva hit oe Gon os lio EC a Se Rr i 1 
ee ale es tec Ie ier elie it a 2 
ee ee Seo tk ONIN, Baas. we we 28 
Petersen 2 QoS big Oe wa er Gio Ss ong nc > Ents ecu nea oa 


SUMMARY 


The relativistic K-shell photoeffect is analysed in terms of four “reaction amplitudes”’, 
which describe the process completely, including all polarization effects. Expansions of these 
amplitudes in radial matrix elements are given. Using an improvement of the Sommerfeld- 
Maue final state wave function, the amplitudes are computed to relative order «Z for a pure 
Coulomb field. Among effects obtained from the first-order terms are a non-zero forward 
emission intensity, and a transverse polarization of the photoelectrons, similar to that ob- 
tained in relativistic electron scattering. 


Introduction 


For a review of theoretical work in the atomic photoeffect up to about 1956 
—which in this context is practically equivalent to 1936— we refer to the 
article by Bethe and Salpeter in Handbuch der Physik XXXV (1957). 

The theoretical interest in the photoeffect has been revived during the last 
years; this is to a large extent due to the more accurate experimental data on 
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angular distribution of photoelectrons which have become available in the 50’s, 
above all from the work of Hedgran, Hultberg and collaborators (see Hultberg 
1959 and references there), who study the photoeffect as a means for indirect 
gamma-ray-spectroscopy, e.g. in the determination of internal conversion coeffi- 
cients. New techniques for the study of electron and photon polarizations, partly 
inspired by the detection of parity violation in weak interactions, have made 
it of increased interest to investigate theoretically the polarization effects in 
the “classical”? electron and photon processes. However, as far as the present 
author is aware, no experimental studies of electron polarization in the photo- 
effect have yet been reported. 

The standard relativistic formula, valid in lowest order in «Z in a Coulomb 
field, for the angular distribution of K-shell photoelectrons from linearly po- 
larized photons was derived by Sauter 1931. The correction of relative order 
aZ was recently obtained by Gavrila 1959. Electron polarizations have been 
studied theoretically in a number of papers: McVoy 1957, Olsen 1958, Banerjee 
1958, 1959, and Fano, McVoy and Albers 1959. These authors in general work 
in the lowest order approximation, the same approximation that gives Sauter’s 
expression for the differential cross section; however, Banerjee also obtained 
those terms of relative order «Z that survive in the high-energy limit. Con- 
cerning the total cross section we can mention—apart from the « Z-correction 
to Sauter’s total cross section obtained in Gavrila’s work—that the complete 
«Z-dependence in the high-energy limit has been worked out, mainly in papers 
by Erber 1959 and Pratt 1960. 

In the present paper we try to analyse the K-shell photoeffect, including 
all polarization effects, in a way which can serve as a convenient starting 
point both for approximate analytic treatment and ‘“‘exact’’ numerical calcula- 
tions. Simple invariance arguments allow one to express the relevant matrix 
elements in four scalar “reaction amplitudes’’, which contain all information 
of the process; these amplitudes then seem to be a reasonable object for study. 

The crucial point in the relativistic photoeffect is the final state wave func- 
tion, describing the outgoing photoelectron; in the Dirac case this wave function 
can be obtained in exact form only as an infinite series of partial waves. With 
a view to numerical calculations the resulting expansions of the reaction am- 
plitudes in radial electric and magnetic multipole moments are given. We then 
also derive expressions for these amplitudes to relative order «Z% in a pure 
Coulomb field. The real parts of the «Z-terms give the correction of order «Z 
to Sauter’s angular distribution, which was calculated by Gavrila. The imaginary 
a Z-terms give a new photoelectron polarization effect, a transverse polarization 
perpendicular to the scattering plane, obtained also when unpolarized photons 
are used, The imaginary terms also contribute to the forward and backward 
emission intensities in lowest non-vanishing order. We also get the « Z-correc- 
tions to the electron polarization effects obtained when the photons are circu- 
larly polarized. 

Gavrila used in his work a final state wave function obtained by a second- 
order Born approximation. We proceed in a somewhat different way and use 
an improvement of the so-called Sommerfeld-Maue wave function; in this way 


we avoid the divergences that occur in second order when one applies the Born 
approximation to the Coulomb potential 
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1. The expression for the differential cross section, including polarizations 


The formula for the differential photoelectric cross section is generally derived 
by first-order time-dependent perturbation theory for the transition from an 
initial state, consisting of a photon with definite momentum and _ polarization 
and an electron in a bound state, to a final state in a continuum, consisting 
of a free electron with a definite asymptotic momentum. The electron states 
are taken to be one-electron states in the central field of the nucleus; the 
influence of the other electrons in the atom can be taken into account ap- 
proximately by the use of a screened potential. The photon is either treated 
as a quantized field or, more elementary, as an external perturbation, harmonic 
in time, in the Hamiltonian of the electron. By standard procedure one then 
obtains the following expression for the differential cross section (= transition 
probability per unit time and unit solid angle, corresponding to an incident 
photon flux of one photon per unit time and unit area), assuming the electron 
to be described by the Dirac equation: 

3 
\. (1.1) 


: eer bay alee tp Cfme\® 
ai fi)-5-(*) he hy 2 <) 


Here p and H=(m'c'+c*p*)! are momentum and relativistic energy of the 
outgoing electron, k, @ and hy momentum, polarization vector and energy of 
the photon. m and e are electron mass and charge, and «@ stands for the Dirac 
matrices. y; is the wave function of the bound electron, normalized to one, 


[oree exp (¢k-r) p,;dr 


fyfydr=1. (1.2) 


w; is the final state wave function, an incoming scattering state with asymptotic 
form 
(exp (tp:r)+inc. spher. wave) w. (1.3) 


-w is a constant 4-spinor, with w*w=1. For simplicity we suppress the loga- 
rithmic term in the asymptotic plane wave exponent that occurs in the case 
of a Coulomb field. 

That an incoming (and not an outgoing) scattering state must be used is 
due to the fact that only a superposition of such states gives for large positive 
times a plane wave packet (without spherical part), describing the photoelectron 
that moves away from the atom in a definite direction. (See Bethe and Sal- 
peter 1957, p. 382, and references there.) 

The factor (mc/h)® in the brackets in (1.1) is inserted to make this part of 
the expression dimensionless. From now on we shall put m=h=c=1, ie. we 
use i/mc, me and mc* as units of length, momentum and energy. 

(1.1) refers to the transition between definite initial and final spin states of 
the electron under absorption of a photon with a definite polarization. To obtain 
the observed cross section we must perform a suitable averaging and summing 
over initial and final state polarizations. We shall give the cross section for 
a general initial photon polarization and a general polarization sensitivity of 
the detector of the photoelectrons. We will not assume the polarization of the 
bound electron to be observed, though this should in principle be possible, e.g. 
by studying the polarization of the X-ray emitted when the hole in the shell 
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B,n, 


Gedy oi 


k,n; 


Fig. 1. 


is filled up; it is, however, not difficult to generalize the intensity formula to 
include also this case. Polarization states of electrons and photons are con- 
veniently described by 2x2 density matrices and parameters related to these 
matrices. We refer to Tolhoek 1956, whose notations we shall follow, for a good 
review of this formalism, Before applying this technique to our case, we must 
study the spin and photon polarization dependence of the matrix element. 

We assume the initial state to be a Sj-state, e.g. a K-shell state. Then both 
the initial and final states of the electron are specified by giving the large 
component 2-spinors u; and u,. (See Tolhoek, loc. cit., p. 279, for a discussion 
of this way of describing the spin state of a relativistic electron.) For the final 
state, uy, means the constant 2-spinor of the asymptotic plane wave part. We 
may then write, observing that y; and y, are linear in u; and u;, respectively, 


J yf a-@ exp (ik-r) ydr=uf M (8) w, (1.4) 


where M(é) is a linear operator in spin space, and can thus be written as a 
linear combination of the unit matrix and Pauli spin matrices. M (é) depends 
linearly on @. 


Introduce unit vectors fi; and fi, in the photon and electron directions (see 


. . wa ~ ~ n x n 
Fig. 1). @ is the angle between n, and ny. hai eA is the normal to the scat- 
i My 


tering plane. As basic photon polarization directions we choose @, in the scat- 
tering plane and @,= 1, so that 6,xé,= fj. 
The structure of M(é) is obtained by observing that if the electron moves 


in a central potential the matrix element (1.4) must be invariant under reflec- 
tion in the scattering plane, 


ugr M (8p) Uin= uj M (8) uw. (1.5) 


Here @, is obtained from @ by changing the sign of the component along @,, 
and uz;=o'nw, apart from an irrelevant phase-factor. 
Then we get from (1.5) 


M,=M (@,)=A (6) +i B(0)o° Ni, 


M,=M (@,) =i 0 (6) 6: (Ax fy) +i D(0) 0: ty. sor 
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A, B, C and D, the “reaction amplitudes’’, are complex-valued scalar functions. 
The factors i are inserted for later convenience to give the same phase to the 
functions in the lowest order expressions for a Coulomb field. 
If 6=0, i. t,=Tfiy (then we define the scattering plane by i=@,), M, is 
obtained from M, by a rotation x/2 around fj. This gives 


A(0)=D(0)=0; B(0)=—C(0). (1.7) 
Similarly A (x)=D(x)=0; B(x)= C(z). (1.8) 


Let us now return to the derivation of the intensity formula. Put é=a, é@, + 
+d 6, where |a,|*+|a,|7=1, and assume the detector of the photoelectrons 
to respond to the electron spin state according to @ae,=4(1+-6). Averaging 
over initial electron spin states (factor 4 in (1.9)) we obtain the following ex- 
pression which comes in instead of the brackets in (1.1) in the general formula: 


4 Trace (0a, VM (6) M* (@)). (1.9) 


We now insert M (é)=a,M,+a,M, and the expressions (1.6) into (1.9); instead 
of the resulting bilinear expressions in a, and a, we introduce the parameters 


&,=|a,|?—|a,|*;} &=a,af+a,a7; &=i(a,az—a,a7). (1.10) 


The result of the trace calculation in (1.9), combined with (1.1), can then be 
expressed as 


CAO ME) OSE (Cok), (1.11) 
; TOOLS Pie ee 
pet G 52 (=) hehyme’ a 


and 
2K (6, €)=4(|Al?+|BlP+|C[?+|D[)+Im (4 BY +CD")o- A 
+&,[4(|AP+|B)—|CP—|D/)+Im (4 B*—CD*)o-n] 
+ & [Im (A C* — B D*) €- (nx fy) +_Im (A D* + BC*) C- ty] 
+&,[—Re (A C* — BD*) €- (fix fi) — Re (A D* + BC*) C- fy]. (1.13) 


(1.13) was derived assuming a complete polarization of the photon, i.e. |§|=1, 
but it is easy to see that the formula is also valid for a partial polarization 
|€|<1. Any 2x2 density matrix g can be written 9@=A/,0;+/ 02, where A, 
A,>0, A,y+4,=1, and g, and g, correspond to pure states (complete polariza- 
tion). Speaking more physically: any statistical mixture can be obtained by in- 
coherent superposition of two appropriately chosen pure states. The intensity of 
photoelectrons obtained from this incoherent superposition of photon polarization 
states is obtained by forming the same linear combination of the intensities 
corresponding to pure polarization states, and (1.13) is then valid also in the 
general case; the linear combination of the density matrices amounts to forming 
the same linear combination of the vectors §. 

It should be stressed that (1.11) gives the cross section per electron; multi- 
plication by two gives the cross section for the K-shell. 
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We shall now look at some special cases of (1.11)-(1.13). 

If the incident photons are linearly polarized, the angle between the polariza- 
tion plane and the scattering plane being ¢, and the final electron polarization 
is not observed (the detector accepts all electron polarizations), we should put 
é=cos Jé,—sin $6, ie. &=cos 24, & = —sin 2¢, and &,=0, and add K (6, &) 
and K(—{€, &). This gives 


o (6, )=C' [C2 +| DP? + cos? g [| A? +|BP—|CP—-|D/I}- (1.14) 
Unpolarized photons correspond to §=0: 
o(0)=C' ¥ {JAP +|BP + CP + DPS. (1.15) 


The same differential cross section is obtained also for circularly polarized pho- 
tons, ¢,=&,=0. 

The polarization P of the outgoing photoelectrons, related to the density 
matrix by 9=4(1+P-o), is obtained from the relation 


pees teres) 
K(t,€)+K(-f, €) se 


We will give P only for the case which is probably of greatest experimental 
interest, namely when the photons are (partially) circularly polarized, €,=&,=0, 
&,=P,. P,,=1(-—1) corresponds to photon helicity +1(—1), i.e. right (left) hand 
circular polarization, P,=0 to unpolarized photons. From (1.13) and (1.16) we get 


PF 


P=P,ii+P,(P, iy + P, t,x hi), (1.17) 
2 Im (A B* + CD*) 
witl P,= 
1 TAPE +[BP+[CP +|DP ee 
pi. al >Re ADpEBg® 
TAP BESTOP +1 DP oe 
p.a_2 Re (AC*-BD*) 
[AF + [BP +P +1 DF re 


Py is a transverse polarization, perpendicular to the scattering plane, and is 
obtained also for unpolarized photons; P, and P, are longitudinal and trans- 


verse polarizations in the scattering plane, and are obtained for circularly polar- 
ized photons. 


2. Initial and final state wave functions 
In this section we shall study the spin structur 1 
' ‘ e of y; 
integral expressions for A, B, C and D. : Yo aoe nee 
y: and y; are solutions of the Dirac equation 


[—ta-grad +a, + V(r) —e]p=0. (2.1) 
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a, %, are the four Dirac matrices, ¢ is the energy (in units mc”). The o-ma- 
trices will be written as direct products of o- and o-matrices (see Dirac 1947, 
p. 255), and correspondingly the 4-spinor space will be considered as a direct 
product of “spin space”, spanned by wu, and w_, and “large and small com- 
ponent space”’, spanned by v, and v_. The o-matrices operate on u, the o-ma- 
‘trices on v. Thus e.g. v_u, stands for the 4-spinor (0, 0, 1, 0). 

The initial S,-state wave function y; has the form 


Yi =[9(r) vx —tfi (r) o- T v_] um, (2.2) 
with normalization 


4nf [git filrdr=1; ue y=. (2.3) 
0 


For the final state wave function we choose the z-axis as the direction of 
p. yy, then has the asymptotic form (cf. (1.3)) 


ue eae 
(F [e+ ine. spher. wave] (v. Se 0.0.) Us. (2.4) 


From the symmetry property of the scattering problem (2.1) and (2.4) under 
rotations around the z-axis and reflection in a plane containing the z-axis it fol- 
lows (see Nagel 1960, Sect. 4) that y, is of the form 


pr {[f+ (7, 0) oF VJs (7, 0) ora Ui + [je (7, 0) th Ye. (7, 0) ora Oz v_} Uf, (2.5) 
where Op = — SIN Yoz+ COS — Oy. (2.6) 


(The use of f and g in (2.2) and (2.5) is somewhat disturbing from the logical 
point of view, but the notation g and f in (2.2) for the large and small radial 
components is standard, and f, and g. in (2.5) are in a natural way connected 
with the scattering amplitudes f(#) and g(%#) in Dirac scattering theory.) 

Having chosen p=pZ, we now put k=k|~—sin 6, 0, cos 6|; then we also have 
é,=|cos 6, 0, sin @| and @é,=n=Yj. 

Inserting (2.2) and (2.5) into (1.4) and comparing with (1.6), we obtain ex- 
pressions for A, B, C and D as integrals over g;, f;, f and g.. The calculations 
can be simplified by using trace techniques and the known forms (1.6). The 
result is 


A(6)=( 1,+1,—I,) cos 6+ (1,—I,;+1,) sin 0 
B(0)=(—1,+1,+1,) cos 0+ (Ip+1,+1,) sin 0 
ObjaeT +1,—I; 
DG) ate Pd — Loe 


Here J, ... J, are the integrals 
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I,= —ifft f, exp (tk-r) cos 0 dr 
I,=—if fifi exp (¢k-r) sin 8 cos y dr 
I,= —if gt f: exp (¢k-r) sin # cos 2pdr 
I,=—ifgth exp (ik-r) cos # cos p dr (2.8) 
I,=—ifg%f, exp (ék-r) sin 0 dr 
l= ffs exp (ik-r)dr 


I,=- foto exp (tk-r) cos y dr. 


3. Expansions of the reaction amplitudes in radial matrix elements 


As is well-known, the exact final state wave function y, cannot, even for a 
pure Coulomb field, be obtained in closed form in the Dirac case. One must 
use a partial waves expansion, and the reaction amplitudes are then obtained 
as series of (ordinary and associated) Legendre polynomials, with certain com- 
binations of radial matrix elements as coefficients. We shall give these expan- 
sions here for future reference, although they will not be used in the following 
sections of this paper. 

It turns out that the expansions in question are simpler if we choose instead 
of B, C and D certain linear combinations, which are related to the circular 
polarization states e, =@,+76@, (photon states of definite helicity) and also relate 
the electron spin to the direction of the photon instead of the outgoing elec- 
tron. Introduce 


M,=M,+iM,; o,=0°0:, o),=6-(8,x6,). (3.1) 


Then we have from (1.6) 


M,=A(6)+E (0) 0, +F (6)0.+@(6)o-, 


M_=A(0)—E(6)o,,—G(0) o, —F (6) o_, a 
where E(6@)=_ sin 0 C(6)—cos 6 D(6), 
F (6) =4[B(6)—cos 6 O(6)—sin 6 D(6)], (3.3) 


G(0)= —4[B(0)+ cos 6 O(6)+sin 6 D(6))}. 


The expansion of yy, fulfilling the asymptotic condition (2.4), is obtained in 
the usual way by taking an expansion, with undetermined coefficients, in partial 
wave solutions obtained after separating the Dirac equation in spherical coor- 
dinates, and determining the coefficients so that after subtraction of the ex- 
pansion of the plane wave part in (2.4) the rest only contains asymptotically 
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incoming waves e'”’, Comparison with (2.5) then gives the following series for 
fz and g.: 


fe (7, 8) = > a [Le g, + (1+ 1) e-"-1-19_;_4] P, (cos 9), 
g(r, “e > # fe g,—e-'"-1-1 g_)_1] Pt (cos 8), 

f-(r, 8) = d [+ 1) e'+1 f.,, Le“ f_1] P; (cos 9), 
g-(r, i > Ue +4 fia —e7'-1 f_] PH (cos 8). 


gn and f, (n=l or —I1—1) are the ordinary large and small component radial 
wave functions, with asymptotic normalization 


e+1\? 1 la 
In © Qe tae pr—-—<+x In 2pr+m, : 


e—1\? 1 la 
h~( ae reat (vr “2-4 1m 2or+ns), 


where 7, are the asymptotic phases, and (if we have a Coulomb potential pre- 
sent) x=aZ/8, «=e/hic, B=v/c. 
Introducing (2.2), (3.4) and the expansion (see e.g. Edmonds 1957, p. 81) 


exp (vk-r) =4x > thi (kr) Yim (8, p) Yin (0, 2) (3.6) 


into the expressions (2.8) for J,...J,, and going over to A, H, F and G@ via 
(2.7) and (3.3), we obtain after somewhat lengthy calculations, which involve 
the use of various relations between Legendre and associated Legendre poly- 
nomials (see e.g. Edmonds 1957, p. 23), the final formulas 


A (0)=42 > {e'™ a, +e'"-!-1a_,_} P} (cos 6), 


H(6)=42 > {—e (+1) [a+ (U1) 6) + e"-!-11 [a_1_1 — (14+ 2) b4- i feos), 
F (0) =22 > {e'™ [(1+ 1) a, + (1-1) by] — "1-1 [La_y_1 + (+ 2) b_1_1]} P; (cos 8), 
G (0) =20 > {fe (0+ 1) [a —b,] +6711 [ — aya t boi ap Flees) (3.7) 


a, and 6; are the following combinations of radial matrix elements, corresponding 
to absorption of electric and magnetic radiation multipoles: 


a, = | ahi 1 (kr) r’ are | tea (kr) rdr— 


5 | hegvies (kr)r*dr; 1l=1,2,... 


=) 
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: 1 . 
I i [a fishies (Er) dr — [ir g-1-1jisi (kr) dr — 


2(1-+1) Ve ae 

— —__— 121 (EDT OTs) t= pee 
yea | feoraieal ) 

b; = | wh thadirken ran [=2, 3, ... 

b_1-1= [ tostratharadina (er) Pdr l=0, 1, ... (3.8) 


a; and a_;_; correspond to electric multipoles of order J. b; and 6-1-1 correspond 
to magnetic multipoles of order /—1 and J+], respectively. 
For the total cross section we obtain from (1.15), (3.3) and (3.7) 


o=0'4 | {14P+|BP+1CP+1D 1 dQ= 


=320°C’ > {(l+1) a? +lary_1 tl bia + (I+ 1) 62}. (3.9) 
l=1 


The quantities a, and 6, that occur formally in F(6) in (3.7) should be put 
equal to zero. 

The expansions (3.7) can, of course, also be obtained directly from (1.4), in- 
serting for @ exp (ik-r) an expansion in electric and magnetic multipoles, and 
for y, the partial waves expansion. Although one then has to make the spin 
analysis again, this method may have the advantage of leading more directly to 
the form (3.2), and also of giving automatically the special combinations a, and 
b; in (3.7). In the method used here this special choice comes in more indirectly, 
e.g. from the requirement that the total cross section (3.9) should not contain 
any cross products. 


4, An approximate final state wave function in a Coulomb field 
We now assume the electron to move in a Coulomb potential 


Vir)=—a/r, a=aZ, a=e'/ihe, (4.1) 


and will derive expressions for the reaction amplitudes to relative order a. This 
means, as will be seen, that we must have expressions for y; and wp, to order a?. 


For the initial 1 S,-state (K-shell state) we have (see e.g. Bethe and Sal- 
peter 1957, p. 155) 


gi=N yp’ ea j= - Ly his (4.2) 
i i ? ‘ l+y is . 


where y= V1—a? is the corresponding eigenvalue of ¢; 1—y is the binding 
energy of the A-shell. The normalization constant WN, is 


y.= (" +) (2a)?7*3\t gis 
Prarie a et 


[1+0 (a? In a)]. (4.3) 
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In the relevant approximation we can put 


i $ 
(=) =a. (4.4) 


Looking at y,;, we see from the fact that the zeroth order approximation is 
a plane wave, that f, and f_ in (2.5) start as a®, whereas g, and g_ are of 
order a. Inspection of (2.8) shows that, apart from the normalization constant 
W;, all integrals are of order a; I,, which seems to form an exception, really 
also starts as a, as is easily seen by using f_~e'”’, g,~e ” and performing the 
integration. 

To obtain expressions for J,...Z, to relative order a, we must then have 
f. and g, to order a and f_ and g_ to order a. The exponent y—1 in (4.2) 
can be equated to zero everywhere except in J,. 

The desired second-order final state wave function can be obtained by a 
second-order Born approximation; this is the method used by Gavrila 1959 in 
his derivation of the differential cross section to relative order a. Due to the 
infinite range of the Coulomb potential, one must apply the Born approxima- 
tion to the Yukawa potential ¢”/r, and afterwards let 2->0. In this limit 
there appear (in second order) divergences in the matrix elements; however, 
as shown by Gavrila, these divergences cancel in the expression for the cross 
section. 

We shall use a different method, consisting of an improvement of the ap- 
proximate Sommerfeld-Maue wave function. This wave function was used by 
Sommerfeld 1939 for a derivation of the lowest order (Sauter) approximation 
of the differential cross section. The infinite range of the Coulomb potential is 
taken into account exactly, and no infinities occur when we improve the wave 
function to take the potential into account to second order. In addition the lowest 
order expressions for the reaction amplitudes will be valid for all B=v/c, thus 
giving a smooth transition to the non-relativistic results. The Born approxima- 
tion is valid only for a<, being actually an expansion in a/f. 

To derive our approximation for y; we go over to the second-order Dirac 
equation by applying the operator ia-grad—%,+V-—e on the Dirac equation 
(2.1); the result is 

[A+p?—-2eV+V*+i(a-grad V)] yp=0. (4.5) 


We are interested in the incoming scattering solution of (4.5), i.e. the solution 
with asymptotic form (1.3). As every solution of the first-order Dirac equation 
is also a solution of (4.5), it follows from the uniqueness of the scattering 
solutions that a solution of (4.5) and (1.3) is also a solution of (2.1), provided 
the asymptotic plane wave spinor exp (ip-r)w in (1.3) is a solution of the free 
Dirac equation, i.e., with suitable normalization, 


+1\# a , 
w= (° “) (v.+ 2 i0-de) U. (4.6) 


(See Nagel 1960 for a more extensive discussion of this point.) 
The Sommerfeld-Maue wave function (Sommerfeld and Maue 1935; Sommerfeld 
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1939, and Bethe and Maximon 1954) is a solution of (4.5) to first order in V, 
and has the form 


y, =e”? (1 Teer: grad) Fw, (4-7) 


where e'”? F is the solution of 


(A+ p?—2eV) p=0 (4.8) 


with asymptotic form e'”*+ inc. spher. wave. 
If we put p=y,+Y., where p is the exact solution of (4.5), we get from 
(4.5) and (4.7) the following inhomogeneous equation for ys: 


[A+ p?-2eV+V?+i(a-grad V)] yp. = 
} 1 
— eh |v" (1 oe A grad) of Qe (a-grad V) a: grad | Fw, (4.9) 
and boundary condition: incoming spherical wave. 
We get the V® order part of y, by solving (4.9) after dropping all terms in 


V on the left hand side and keeping only V® terms on the right hand side; 
this gives 


me lS (exp (—tpir—r' |) eles ‘ls : ; 
Y2 rel [r—r’| € V Ayt5,(¢ grad'V)(o-grad’F),|dr’w. (4.10) 


F, and (grad Ff’), means that we shall take only the lowest (zeroth and first) 
order term of F and grad F, respectively. 
For V=~—a/r the relevant F is 


eb™D (L+ix) F,(—ix; 1; —ip(r+z)), (4.11) 


where iF, is the confluent hypergeometric function. 

_ Combining (4.7) and (4.10), with F as given by (4.11) and w as in (4.6) 

(p=2Z), and using 

8) eee 
p(r+z) 


(grad ,F,),=—ip(? + (4.12) 


we get by comparison with (2.5) the following expressions for f* ete. to the 
necessary orders in a: 


fLaN feo Flin 1; ip(r+z))—a ae Lap (sy eile 
2(e+1) pr ? 

t=N,|—a PP .-ine EXP [Yp(r+z)]—1 sind 

935) ae east 1) be pion 16 


—ipz 


| RS : : noe : 
fi-v,| Pie ibe FP, (ix; 1; +P) Saat (a2; 1; vp (r+z)) (1+cos #) 


p exp (ip|r—r'|) e- t exp [ip(r'+2')]—1 
Bee | exp [ip 
4a (e+1) |r—r'| y'? TS pr’ dr ,/ 
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g=N,)— 0 oF (tx; 1; ¢p(r+z)) sin 0+ 


Maco cP exp (ip|r—r'|) ¢ ® exp [i p(r’+z')]-1 sin 9 cos (y'—g) , , 
tate 3 . - dr’>. 

8a(e+1) |r—r’| r pr 1+cos 8 
(4.13) 


+ 
Here N= (5) eb™ T (1—ix). (4.14) 


It should be observed that in (4.13) the terms containing the ,F,-function (or 
its derivative) are correct for all x=a/f, whereas the other terms, which give 
contributions of relative order a in the matrix elements, are valid under the 
assumption x<1; only then is (4.10) a good approximation, and it is then 
meaningless to give a more accurate form for the other terms giving contribu- 
tions of relative order a. 

The factor sin # cos (y’—q) in the integral for g* may require a comment: 
from (4.10) one obtains o, sin #, which is rewritten, using 


oy =a: =a: [cos (p —g~) P+sin (g’—G)zx@]. 


The second term gives zero when we integrate over r’. 


5. The reaction amplitudes to relative order a Z 


Inserting (4.2) and (4.13) into the integrals in (2.8) we get a number of rather 
complicated three-dimensional integrals. The evaluation of these integrals, which 
ean all be expressed in elementary functions, is performed in the Appendix. 
We give here the resulting expressions for A, B, C and D in the form 


A (6) =K {A, (0) +a [Air (0) +4Az; (6)] +0 (a®)} ete., (5.1) 


‘where Ay, Aig and Aj; are real functions, and K is a factor common to 4, 
B, C and D. The expressions for A;p, Ay; etc. have been obtained under the 
assumption a<f; this fact has been used when we extract the factor K also 
from these terms. 

For shortness we introduce auxiliary functions X, Y, Z, U and J: 


X=2ek(1—£ cos 0) +2¢(e—1)(1—8 cos 8), [as k=e—-Vi—a’] | 


Niel peor. 008 6, 


e—J 
z—me(ll-B cos §)+ cos 6 In ¢ (1+) 
- (e—1) p? sin® 6 , (5.2) 
id Xt—(e—-1)Y 
~ (e—1) p? sin? 6’ 
Vy 2 Ine(1+ 8) 


openly 
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The results are then 


p sind, 43)-X | 
ST = eer xX? [2 ( ) |; 
3 » sin® 6 
B, = —422(e—1) xX? ? 
(5.3) 
eng fin 0 
Cy) =40 2p (6-1) “as 
sin 6 3 l 
D,=4a re [X—2(e—1)]. 
MT el tl ae 1 
Ain=** sin 6 | 8 xXhieg xt cos Cao ; 
=—1" 1". Ze sl ecosie 1 
= 772 aS, + : 
ey besos p Xt | Xs eal (5.4) 
op wl See ea a 
7 eG (Pe =. 
Cir | PEM Re XxX? “|, | 
Dyp=2' sin 6U. | 
Ai;= —xa sin 0[V+(e—1—p cos 8) Z], | 
Byy= n| (=? cos 6) V+ (2 cos B+ p sin’ 6) 2], | 
p 3 (5.5) 
Cy =n (14+ cos 6s) ¥-@-(e=1) sin? 0) 2], 
Fea eae | 
Dyy=x sin 6|—P_ V+ (e-1) YZ). 
oat: | 
For K we have 
K=1aN, N,(1—ix) et pix in (1- cos 6) exp (- arctg ore} (5.6) 
“n'—€& 


i 
Here 6 is an angle-independent phase-factor, which goes to zero as p—0; the 
value of arctg should be chosen in (0, z). 

The lowest order terms (5.3), which correspond to the Sauter approximation, 
are rather simple, and give reasonably short expressions for the cross sections 
and polarizations. The first-order correction terms (5.4) and (5.5) are much 
more complicated, and we will not give the corresponding analytic expressions 
for the observable quantities. The expression for the first-order correction to — 
the differential cross section, involving products of zeroth and first order real 
terms, has been given by Gavrila 1959. Our relative differential cross section 
should agree with Gavrila’s; we have not checked this directly, but, as will be 
seen later on, the expressions for the integrated cross section are the same. 
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To get an idea of the influence of the first-order terms Aj, etc. on the 
cross sections and polarizations, we have computed numerically the functions 
A, (6) ... Di: (9) for three energies, e=1.5, 2 and 3 (B =0.745, 0.866, and 0.943), 
ie. photoelectron energies of about 0.25, 0.5 and 1 MeV. The results are dis- 
cussed in the next section. 

For small £, Ay is the dominating term; B,, C, and D, are of relative order 
B° at least. Even the first-order terms are of order 6 (or higher) compared to 
A,. As has been pointed out earlier, the expressions for the first-order terms 
are unreliable when 6<a, but an exact investigation, using the expansions in 
Sect. 3, of the limit 60 shows that the corrections of order a actually vanish 
in this limit, so that for f=0 the relativistic corrections are of order a? and 
higher (see Nagel and Olsson 1960). 


6. Results and discussion 


a. Angular distribution, unpolarized photons 


Although the angular distribution (including the case of linearly polarized 
photons) to relative order a has been given by Gavrila, we shall give a dis- 
cussion here, as the formula in Gavrila’s paper is rather complicated. 

The relative angular distribution is (see (1.15)) 


o (0)~I, (0) +al,(6)+---, (6.1) 
where I,=Ao+ Bot 0+ Dj= ess 
1 sin? 0 €(e—1) (e-2) ry 
= 3 aS =* ae » 
32 2 = (e—1) (1—f cos 6)! j1+ 9 (1—f cos @) 


s the Sauter cross section, and 


I, =2 (Ay Ain t+ By Biz t+ Cy Cin + Dy Diz). (6.3) 


In (6.2) we have put k=e—1, cf. (5.2). 

Graphs of J, and J, for e=2 are given in Fig. 2, p. 21. For «=1.5 and 3 
similar curves are obtained, except that the relative importance of J, is smaller 
for «=3 than for 1.5 and 2; see Fig. 7, p. 22, for the ratio of the integrated 
cross sections J, and J). It is seen from Fig. 2 that the influence of J, on the 
relative cross section is very small, both because J, is considerably smaller than 
I,, and because J, has about the same angular dependence as Jj. There is a 
slight tendency of J, to push the angle 6,,, of maximum emission towards 
larger angles and to increase the emission for angles larger than Ox. Both 
these features are present in the experimental angular distributions obtained by 
Hultberg 1959, but especially the emission for large angles, obtained in the ex- 
periments, is much larger than can be explained even qualitatively by the first- 
order correction J,. For heavy atoms—Hultberg’s results refer to Z=92—the 
higher order corrections in a@ must be expected to be at least as important as 
the first-order correction. 
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From (6.3) and the fact (see (5.3)) that Ao, Bo, Co and D, vanish for 6=0 
and a, it follows that the forward and backward intensities must be zero even 
in first order (cf. Gavrila 1960). The forward (and backward) intensity is of 
order a?, and the a? order term is evidently obtained in the form (cf. (1.7)) 


I (0) =a? 2 (Bip + Bi). (6.4) 
In Fig. 3, p. 21, we show the resulting ratio 
%y=1(0)/Ip (Omax) (in units a”) (6.5) 


of forward to maximum emission, as function of energy. The figure also shows 
the three experimental points obtained by Hultberg for Z=92. The theoretical 
and experimental ratios show the same tendency to increase with energy, but 
the quantitative agreement is bad; again, this should not be surprising, in view 
of the importance of higher order terms for such large Z-values. Banerjee’s 
conjecture (Banerjee 1958) that x) should decrease as the energy increases was 
based on the incorrect assumption that the finite forward intensity should be of order 
a/e. Sauter and Wiister 1955 have performed a numerical calculation of x9 for 
hy=1332 keV and Z=82, obtaining the value 0.27 in rather good agreement 
with the experimental 0.245 (for Z=92). They also give an approximate value 
0.21; this value was computed taking into account only a part of the first- 
order corrections in the matrix elements, which explains the disagreement with 
our result (~ 0.42). 

We also obtain in order a? a non-zero value for I (z)/I, (Omax), but this ratio 
is very small, 0.007 a® for ¢=1.5, and decreasing rapidly as the energy in- 
creases. 

It is possible that the vanishing in lowest order of B and C for 6=0 and 
za is accidental for the Coulomb potential. If we use a screened potential we 
might then get a finite forward intensity even in lowest order. This could be 
so only for relativistic energies, because non-relativistically, when spin effects 
are ignored, B, C and D must vanish identically (cf. (1.6)). 


b. Angular distribution, linearly polarized photons 


This distribution is obtained from (1.14). If we take the lowest order (Sauter) 
approximation, for ¢=1 (low-energy limit) no emission takes place perpendicular 
to the polarization plane, but as the energy increases, part of the electrons are: 
emitted also in the perpendicular plane, until for ¢>1 the distribution is in- 
dependent of the azimuthal angle ¢. The effect of the terms of first order in 
a is rather small; in the energy range ¢e=3 to 5 they increase somewhat, by 
a factor 1+0.3a, the ratio of emission in and perpendicular to the polarization 
plane. The behaviour for ¢=1 and for e=oco is not changed. 

If one studies the asymmetry ratio for a fixed 6, say z/2, there even occurs, 
in the Sauter approximation, a ‘‘crossover’’, so that for e>2 the emission in 
the polarization plane is smaller than in the perpendicular plane (see graph in 
McMaster and Hereford 1954), The effect of the a-order correction is to di- 
minish this ratio, by a factor which is about 1—a in the region e=1.5 to 3 
thus causing the crossover to occur at a somewhat lower energy than e=2 
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The experimental results on this point are somewhat contradictory. McMaster 
and Hereford, loc. cit., using Compton-scattered photons to produce photoelectrons 
from Pb, found the crossover to occur in fair agreement with the Sauter re- 
sults, whereas Brini et al. 1957, using an experimental arrangement, which should 
enable them to count only the K-shell photoelectrons from Pb, found the ratio 
to be larger than one everywhere in the interesting energy region. 


c. Electron polarization 


We shall study only the polarizations P,, P, and P3, (1.17)-(1.20), obtained 
from circularly polarized or unpolarized photons. Let us first remark that from 
the relations (1.7) and (1.8) it follows, if B(0) and B (z) are different from zero, 
that P,=P,=0 for 6=0 and a, and P,(0)=1, P,(a2)=—1. However, in the 
approximate expressions studied below, these relations may not be fulfilled, as 
a consequence of the vanishing of B, (and C,) in the forward and backward 
directions. 

In the limit #0 all polarization effects vanish in the approximation con- 
sidered here. However, the higher orders in a give non-vanishing and for heavy 
atoms even considerable polarization effects also at the K-shell threshold, as 
will be seen in a following paper (Nagel and Olsson 1960). 

P,, the transverse polarization perpendicular to the scattering plane, is a first- 
order effect in a, and is given in this order by 


oA eB er Al; Bi Cs Dis + Or DI 


a rh Ag+ Bh + 03+ Di 


(6.6) 
Graphs of P,/a as function of @ for the three energies are given in Fig. 4, 
p. 21. P, is, for these energies, fairly constant over the angular region where 
most of the photoelectrons are emitted. P, is a polarization analogous to that 

obtained in relativistic Coulomb scattering of unpolarized electrons, and can be 
studied experimentally in the same way by observing the left-right asymmetry 
in Coulomb scattering of the outgoing electrons. The transverse polarization in 
electron scattering vanishes as ¢' for large e; as will be seen in (e) in this 
section, P, also behaves, at least in first order in a for the photoelectrons going 
out in the peak of maximum emission, as e', so that it must attain a maxi- 
mum for some energy larger than ¢=2. The higher orders in a give a P, dif- 
ferent from zero even in the limit 6=0. 

The longitudinal and transverse polarizations P, and P, in the scattering 
plane, obtained with circularly polarized photons, are different from zero even 
in the lowest order approximation: 


— p 1+the(e?—2e—1)(1—B cos 6) 
2 e+11+4e(e—1) (e—2) (1—f cos 8) 


sin 6 
1+2e(e—1) (e—2) (1—f cos 6) 


P,=(e-—1) 


In Figs. 5 and 6, p. 22, we give graphs of P, and P35, as obtained from 
(6.7) and (6.8), and also of the first-order corrections { (8), defined by 
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Proorr = [1 +f (9)] P. (6.9) 


For the longitudinal polarization P, the effect of f(@) is to diminish the po- 
larization by about 1—0.3a in the angular region of greatest interest, for these 
energies. 

The correct lowest order expressions for the electron polarizations corresponding 
to circular photon polarization were first given by Olsen 1958. Banerjee (1958, 
1959), using the Sommerfeld-Maue wave function, arrived independently at these 
results, and also gave those corrections of order a that survive in the limit 
eco. In a paper by Fano, McVoy and Albers 1959 there is a very extensive 
and illuminating discussion, with many graphs, of these lowest order polariza- 
tions. 

To our knowledge no experimental investigations on electron polarizations in 
the photoeffect have been published. 


d. Total cross section 


The total cross section (per K-electron) is obtained by integration of (1.15), 
and can be written in the form 


o=Tos(l+af(e)+-:-), (6.10) 
. 2 (1+x?) #(e+1) 
with T = Loe &XP | -2% arctg micas | cf. (5.6). (6.11) 
ds is the Sauter cross section 
8x (nh\? mc?\* 
Oa Tae (=) a® Z* (7) p Ny (e). (6.12) 


_f4,0)do Y,(e) 


f (e) si, (4a Male) (6.13) 
with 
= 3 (e—2) tal 1+ 8 
No le) a) eaeneaa (1 spa 75h fees 
3a pp —192*+41le+8 2 (¢—2) 1 l= 
N 
1 (€) v1 | i (e=1) eerbakesail (6.15) 


a the ee of cae we have used, in the dominant term for small B 
e exact expression for X in (5.2); this gi f i 
een a (5.2) 8 gives four of the five powers of | 
The formula for o agrees with that obtained b i i 

y Gavrila 1959, if we expand 
T to first order in x. For small 8, 7’ os goes over into the non-relativistic Stobbe 
cross section (see e.g. Bethe and Salpeter 1957, p. 390). f (e) goes to zero with f; 
this is strictly true, even if the form given above for f(g) is not reliable when 
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Bsa. A graph of f(e) as function of ¢ is given in Fig. 7, p. 22. The behaviour 
of f(e) for small f, and for large e is: 


B<1, f(e)= 5B (1+ 6*/2+--), (6.16) 
e>l, He= -FE(1-F ++), (6.17) 


For large e, f(e) goes to —42/15= —0.838. This limit is, however, approached 
very slowly, as is seen from (6.17) and Fig. 7. 

To give an idea of the accuracy of (6.10) we have tabulated (Table 1) the 
ratio o/os for three values of Z and three photon energies (Z=26, 50 and 84; 
hy/mc’ =0.693, 2.212 and co), compared with the “exact” values in ( ); for 
the lower energies these exact values are taken from the numerical calculations 
of Hulme et al. 1935, for hy/mc*=co from the work of Pratt 1960. The values 
in [ ] are obtained from an approximate formula proposed by Pratt, which 
combines the a-dependence of his high-energy limit cross section with the energy 
dependence given by f(e) (from Gavrila’s work). This formula should be a good 
approximation at high energies. 


Table 1 

hvime | Z=26 | Z=50 | Z=84 
0.59 0.45 0.44 

0.693 (0.62) (0.50) (0.51) 
[0.60] [0.44] [0.34] 

0.60 0.43 0.32 

2.212, (0.63) (0.51) (0.47) 
[0.64] [0.50] [0.44] 

ies 0.49 0.26 0.11 
(0.53) (0.34) (0.22) 


e. The high-energy limit 


For large ¢ the emission of photoelectrons is concentrated in a narrow peak 
around an angle 6§ of order e¢*. We will keep only the dominant terms in ¢ 
of A, ... Di; for these 6, and study the resulting expressions, However, it should 
be stressed, cf. the behaviour of f(e) in (d), that these limits can be expected 
to be approached very slowly. 

Put x=ec6; then we have e.g. X=1+<2" ete. After suitable expansions we 
get to first order in a: 

A;?+|B/?=|CP?+|DP=|KP (8 pact beter (6.18) 

ais oe. whaed 0) coy ir of 2Y bead]? 
This shows, if we compare with (1.14), that to first order in a the distribution 
is, for these 9-values at least, independent of ¢ in the high-energy limit. 
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(6.18) attains its maximum 


; 4 I 


1 It 
2=— 1+a—7-). (6.20) | 
me “al 8V6 : 
Thus, even in the high-energy limit we still have the tendency of the first- : 
order correction to increase the angle Onax, though the correction now subtracts 5 


from the Sauter distribution. 
For the forward to maximum intensity we find the asymptotic value 


Xo = 5.85 a*. (6.21) | 
For the polarizations we obtain, to order a, 
P,=0(a/e), P,=1+0(e), P3=0(e7’). (6.22) | 


The high-energy results given here agree, in the relevant cases, with those that 
can be derived from the formulas given by Banerjee 1959. 


7. Conclusions 


Some of the effects studied in this paper, as the finite forward emission and 
the transverse polarization of the photoelectrons, come from the terms of rela- 
tive order «Z in the matrix elements. For the other effects, which appear 
already in the lowest order approximation, the corrections obtained here are 
in general rather small, at least for Z-values (Z2<30-— 40, say), where one can 
hope that the first-order terms give the main part of the corrections. For heavy 
atoms, on the other hand, the first-order terms can at best be expected to give 
a qualitative indication of the effect of the corrections. To get reliable results 
for the higher Z-values one has to use the expansions in Sect. 3. This means 
that one must compute numerically a large number of radial matrix elements, 
especially if one is interested in the differential cross section and polarizations; 
for the total cross section relatively fewer radial matrix elements will be ex- 
pected to suffice. This numerical computation can probably be carried through 
with a reasonable amount of work only in the case of a pure Coulomb field; 
in this case, the formalism and series expansions worked out by Olsson 1959 
in a number of papers should be of great value. At the K-shell threshold, where 
the outgoing photoelectrons have zero kinetic energy, one should expect the 
series expansions in radial matrix elements to converge rapidly, even for the 
heaviest atoms. This limiting case is studied in a following paper (Nagel and 
Olsson 1960). 

Even if these exact calculations are carried through, it must be remembered 
that they assume the photoelectron to move in the field of a bare point-nucleus. 
Although it is probable that, for not too low energies and Z-values, the in- 
fluence of screening on the total cross section and the general form of the 
angular photoelectron distribution is not very large, such effects as certain de- 
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ails of the angular distribution—e.g. the forward intensity—and the polariza- 
ions may be more sensitive. Other corrections, as the influence of the finite 


uclear size and radiative corrections, are probably much less important than 
he screening corrections. 


ACKNOWLEDGEMENTS 


In the first place I want to thank Professor L. Hulthén for facilitating my work in various 
ays. Dr. S. Hultberg’s experimental investigations formed one of the main inspiration sources 
or this work, and Dr. P. O. M. Olsson has contributed much, both in discussions and by his 
nthusiasm. My thanks are also due to Dr. U. Fano and Dr. R. H. Pratt for kindly sending 
e preprints of their papers. 


The financial support from Statens Rad for Atomforskning is acknowledged with gratitude. 


Royal Institute of Technology, Stockholm, Division of Mathematical Physics. 


2000 - 
tops 
3 
2 
1 
t 
Sa a PS 8! a -— 
30 60 90 co) 1 2 3 L 5 a: 
Fig. 2. Fig. 3. 


Fig. 2. The lowest order (Sauter) angular distribution J, (@) and the first-order correction J, (0), 
for €=2 (hv=0.5 MeV). See (6.1). 

Fig. 3. The ratio x, of forward to maximum intensity, in units a*, as function of energy. See 

(6.5). The three experimental points, for Z = 92, were obtained by Hultberg 1959. 


F/a 


Fig. 4. The transverse electron polarization P, (in units a), perpendicular to the scattering plane, 
as function of 0, for €=1.5, 2, and 3. See (6.6). 
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Fig. 5. Fig. 6. 


Fig. 5. The longitudinal electron polarization P,, corresponding to circularly polarized photons 
in lowest order, and the first-order correction f (0), for three energies. See (6.7) and (6 9). 
Fig. 6. The transverse electron polarization P, in the scattering plane, corresponding to circu- 
larly polarized photons, in lowest order, and the first-order correction f ((), for three energies. 
See (6.8) and (6.9). 


Fig. 7. The ratio f(¢), (6.13), between first and zeroth order total cross sections, as function of 
energy. 
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APPENDIX 


In this appendix we shall indicate the computation of the integrals J, ... I,, 
obtained by inserting (4.2) and (4.13) into (2.8). 

To save space we use the following abbreviations in (A. 1) and (A. 2): 
F=,F, (tx; 1; 1p(r+z)), H=exp [—art+i(k-r—pz)],J= texp [i p (r+ 2)]— 1}/pr. 
It is pomverient to write every integral as a sum of two parts, J° and J’, as 

~ follows: 


N=Ni| HF 008 dx | 
dal sin #? cos pdr 


w= 


“| arar—[ HP 1+ cos #)d | 
ei 


Rane | WP sin # cos pdr 
u 


1 ae 
re. at hoe d 
fi Aceon pi | 43 cos ddr 


SNe aceng i [ HI sin 9 cos par 
oe chal sin’ # cos 29 
a= Creal Theos et 
ae pa’ pein @ cos O cos @ 
He Nas Pi | ad 1+ cos 8 
ee, 
oats p ‘ sin” ? d (A. 2) 
sie Nas Poi | nase cos } ; 
ee P 4 | a-P)H Pdr 
Ig= Na sats ( 7 ) 
aad "ly Hir)¢?¢-*? ft(tes ~J(r r))drar'} 
m, |r—r’| 
expip|r—r'| 5, ip (z—2’) ,./-2 Aye 
7 re r’-* J (r’) 
a Shi pian senas || 7 fr—r'f His 
sin # cos @’ ; 
. -ardr’, 
1+ cos 3” J | 
Here N=iaN,Ny, (A. 3) 


where N, and N; are given by (4.3) and (4.14). 
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As the integrals J’ give contributions only in the first (relative) order in a, 
we can, at least in the final results, put a=0 in the terms H that occur in 
(A.2). The factor cos g’ in the expression for I} is obtained after some sim- 
plification: first one obtains cos (p—g’) cos 9; the corresponding integral with 
sin (p—q’) sin g evidently has the value zero, as the integrand is an odd 
function of the variable g—g’, so we can add this integral, and hence obtain 
cos g’. 

The integrals 7°, which give the lowest order (Sauter) expressions for the 
matrix elements, can easily be computed from the integral 


I (a, q, p=| exp (—art+iq-r)r',F, (iu; 1; i (pr+p-r))dr= 
(A. 4) 


4a ( a+¢ i 
eg \att+gt+2p-q—i2pa) — 


Here (for a real and >0, p and q real, and p>0) one should use arg (a* + q”) =0, 
and arg (a®7+q7+2p:q—i2pa) in (—7, 0). 
This integral is given by Sommerfeld 1939, p. 462. It is a special case of 
an integral computed by Nordsieck 1954; see also Nagel 1960, Appendix A. 
With J (a, q, p) as defined in (A. 4), and 


qg=k—p=|-—k sin 6, 0, k cos 0—p|, p=pzZ, (A. 5) 


we can write (cf. Sommerfeld, loc. cit., p. 487) 


0 A C mz 
R=N|Sraq+tt,y)| , B-x[Sr@q+tep] 
ot t=0 ot t=0 
1 Ze @ ) 
P= N= — ————_- — res A 
6 “| e+ bat q; Distt a Aas a. va}, (A. 6) 
ay aaa 
7 |Sre q, PZ +3] 


Performing the operations indicated in (A. 6), we get for the parts of A, B, 
C and D corresponding to I°, including first order terms in a: 


A® (6) = B sin 6 [e +3 —e (e—1) (1—B cos 6) +0 (a?)] 
B® (6) =B[—p (e—1) sin? 6—iae (e—1) (1—f cos 6) +0 (a’)] 
C° (0)=B[p? sin? 6+iae(e—1) cos 6 (1—B cos 6) +0 (a?)] (aeg 
D° (0) =B sin 6 [p (e—1) — p? cos 6+iae(e—1)(1—B cos 6)+0 (a?)], 
2 2 ix 
where B=8a2(1—-ix w( tea : P 
N\G iG +2 pq: 2 —i2pa @rgerr 
For the expressions in (A. 8) we have 
X=a>+q=2ke(1—f cos 8), (A.9) 
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a+q+2pq-%Z—i2pa eb ix “(l ee (24 10) 


(A. 10) gives via (A. 8), and ignoring order a”, the three last factors in (5.6). 
The computation of the integrals Jj ... I; in (A. 2) is more complicated. We 
start with J} and J}. 
Integrating 


[exp (-wrtiqrirtar=dayesey (Re «> 0) (A. 11) 


with respect to uw, we obtain 


: é 47 27% tut 
exp [—urtiq:r]r*dr=— arct = In 
{ p[-prtiq-t] , gq/u oo aes 


(A. 12) 


Applying twice, with ~=a and w~=a-—ip, the trick used to derive I} and I? 
in (A. 6), we have 


: 2a Vitek 
lim err keep) celeste) e-rar= 
a0 pr =m es 
moet l <9 ~ i ptk 2pk).a a : 
a Pat Beata jn ea ear e-kr, q=k—pz 

(A. 13) gives 

D=Na i kee in R cos 6]. (A. 14) 

: e+1 X! 

h=Na sin 6 |? bya tiak], (A. 15) 
— In e(1+) 
where X and Y are given in (5.2), and Re" a a B : 


Postponing 73 till later, we turn to Ti and If. In Jj we decompose 


sin # cos B cos p _ sap lose _sin # cos ~ 
1+ cos # ?~ "T+ cos 3 | 


The first term gives J3, while the integral over the second term is equal to 


‘n sin 6 
lim 2 Etie: pt +19)| er vi me TE (A. 16) 
x30 2% Ot t=0 
as is seen from the relation 
rin ‘ exp [ip(r+z)]—1, ,_ 2 J 
Fi (ix; 1; ip(r+z2))=% prea — +0 (x"*). (A. 17) 
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Then we get 
fe : g5 1 -22 5)| A. 18 
L=Na sin 6 | Pers qitia(R eT x ( ) 
nO 


To compute J} we observe _8in UV =1—cos%. The second term is J}, while 
1+ cos # 


the first term is seen to give 


Pe Z)=- — A. 19 

ae I (a, q, pZ) 4upy- (A. 19) 

pl = pe—1 Y ; 2p ¥ A.2 

Hence I va| Ts t gat ize (R cos pice ix)l: (A. 20) 


We shall now compute Ij, and start from 
sin? 9 cos 2y= sin? #—2 (sin @ sin ¢)’. 


The first term gives Jj; for the integral over the second term we obtain (ob- 
serve that q: ¥ =0) 


. } afl A 
[coca pe 


o 


~~ | aafiim? | 1a +sj, patty = 
x0 # cs ot a YP y) s=t=0 P? 


BS (A. 21) 


Altogether, we obtain for J} 
hed? cu ; 2 ne 
Then ltt 205 Ge aA Py ty 
3 als ( eid yt a) tia(R cos hee 2Z}|. (A. 22) ¢ 


U and Z are defined in (5.2). 
For the first integral in Jj we have in the desired order 


[a-rynrarayet | fm rH (r)dr| re (A. 23) 


The integration over angles is easily performed, and in the integration over r 
one can rotate the path of integration to the imaginary axis: 
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47 2 a0" 
| In rH ede" Im | rin r- exp | —(a@—29)rlar= a0 (a, (A. 24) 
q 
0 


The first integral in 75 then gives a contribution 


2p 1 
Nag ee 
ant — = (A. 25) 
I; and the second integral in I} are of the form 
n= | [ex —ar+ik-r) eee eae a f(r’)drdr’. (A. 26) 


From the Green’s function property of exp (ip|r—r’|)/|r—r’| it follows that 
I (f)= ~4x{ exp (—ar+7k-r)u(r)dr, 


where u(r) is the solution of (A+ p*)w=f, which for large r behaves as an 
outgoing spherical wave. Using Green’s formula, with w as above and v= 
=exp (—ar+ik-r), observing that (A+k*’)v=0, apart from terms of order a, 
we obtain 


Li= -e [ ro exp (—ar+ik-r)dr+0 (a). (A. 27) 


This result could also have been established by the use of a suitable Fourier 
representation of exp (ip|r—r’|)/|r—r’ |. 

The integrals in J; and J} can now be computed in a way analogous to the 
earlier computations; no essentially new difficulties arise. The results are 


‘a Pap ont | ee =e | ine 
h=Nalx (2. 5 part 2 +14 1) pea A B cos @)|7, (A. 28) 


In e (14+ 8) 
2 


H=Na sin 6| —2°U +i —(s—1) ra)]. (A. 29) 


From the expressions for Jj ... J; and the a-order terms in A®, B®, C® and D® 
we now obtain Aj;z ... Di; in (5.4) and (5.5). 
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